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Abstract. In this paper we start the study of Schur analysis in the quaternionic 
setting using the theory of shea hyperholomorphie functions. The novelty of our ap- 
proach is that slice hyperholomorphic functions allows to write realizations in terms 
of a suitable resolvent, the so called S- resolvent operator and to extend several results 
that hold in the complex case to the quaternionic case. We discuss reproducing ker- 
nels, positive definite functions in this setting and we show how they can be obtained 
in our setting using the extension operator and the slice regular product. We define 
Schur multipliers, and find their co-isometric realization in terms of the associated de 
Branges-Rovnyak space. 



In this paper we develop Schur analysis, and in particular the Schur algorithm, and a 
theory of linear systems when the complex numbers are replaced by the skew-field of 
quaternions. An important tool is the theory of slice hyperholomorphic functions. So 
there is a combination of a non- commutative setting (since the quaternions lack the 
commutativity property) and of analyticity (via the slice hyperholomorphic functions). 
Since the paper is aimed at two different audiences, namely researchers from Clifford 
analysis and researchers from operator theory and classical linear system theory, we 
will survey the basic definitions from both fields needed in the paper. 
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We denote by S the set of functions analytic and contractive in the open unit disk 
D C C. Such functions bear various names, and we will call them Schur functions 
in the present paper. Let s G iS, and assume that |s(0)| < 1 (and therefore, by the 
maximum modulus principle, s is not equal to a unitary constant, but takes strictly 
contractive values in D). Then, the function 



1 s{z) - s(0) 



1 - 

also belongs to S. More generally, the recursion 

= s{z) 

1 - s(")(0) 



z^O, 
z = 0, 



- |sW(0)|2' 



z = 0. 



defines a sequence, finite or infinite, of Schur functions s^^\s^^\ . . .. The sequence is 
infinite if 

s(")(0)e©, n = 0,l,..., 

while it stops at rank n if |s^"^(0)| = 1. The construction of this sequence is the cele- 
brated Schur algorithm, developed by I. Schur in 1917. See [721 [731 EH] . The numbers 
Pn = s^"^(0) are called the Schur coefficients associated to s, and the sequence (finite 
or infinite) of Schur coefficients uniquely determines s. 

Let us already make the following remark at this point. In the Schur algorithm one 
makes use of Schwarz's lemma and of the elementary fact that if u and v are in the 
open unit disk, so is 

(1.3) 



L — UV 

Let us denote by B the open unit ball of the quaternions. If one replaces in (11. 3p u and 
V by quaternions of norm strictly less than 1, then both 

{u — v){l — uv)"^ and {1 — uTi)^^ {u — v) 

are still in B. But, as we shall see, these transformations will not keep the property of 
being slice hyperholomorphic, and a different approach will be needed. Note that the 
Schwarz lemma will not help to develop a Schur algorithm here because of the lack of 
commutativity. If s is slice hyperholomorphic into B (and in particular |s(0)| < 1), the 
functions 

or 



{s{p) - smil - s{p)s{0))-' 

are indeed contractive in B, but they will not be slice hyperholomorphic. So there is 
no direct counterpart of (II. ip . One needs to use the notion of Schur multipliers and 
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the slice regular product. 

Schur analysis originates with the works of Schur, Herglotz, and others (see [55] for 
reprints of the original works), and can be seen as a collection of questions pertaining 
to Schur functions and their various generalizations to other settings. Among the 
problems we mention in particular: 

(1) Classical interpolation problems such as Caratheodory-Fejer and Nevanlinna-Pick 
problems, and their matrix- valued versions; see for instance [511 [52]. The case of 
boundary interpolation is of special importance. See for instance [26] for related recent 
results. 

(2) Realization of Schur functions in the form s(z) = D + zC{I — zA)^^B, where the 
operator matrix 

'A 

is subject to various metric constraints, namely coisometric, isometric and unitary. See 

[2H112H11S]. 

(3) Schur functions are closely related to the theory of linear systems. The term linear 
system encompasses a wide range of situations. Here we have in mind input output of 
the form 

n 

Vn ^ ^ ^m^n—mi ^ O5 1; • • • ; 
m=0 

where Sq, Si, ... is a sequence of matrices of C^''^^^ (the impulse response), {un)nmo 
a sequence of vectors of (the input sequence), and (|/n)neNo is a sequence of vectors 
of C*"'^ (the output sequence). The function s{z) = Schur function if 

and only if the £2 norm of the output sequence is always less or equal to the £2 norm 
of the input sequence. In other words, Schur functions are the transfer functions of 
time-invariant dissipative linear systems. 

(4) Yet another direction of research is related to inverse scattering; see [321 HH El [7] . 

(5) Last but not least we mention the connection with fast algorithms; see [701 IE] • 
These various questions make sense in more general settings, of which we mention 
in particular the several complex variables case, the indefinite case, the time-varying 
case, the non-commutative case, the case of compact Riemann surfaces, and the case 
of several complex variables, to name a few. References are given in the last section of 
the paper. In the present work we present a counterpart of Schur functions and of the 
Schur algorithm in the quaternionic setting. 

The main tool that we use to extend Schur analysis to the quaternionic setting is the 
theory of slice hyperholomorphic, or slice regular, functions. Some references for this 
theory of functions are [SSI [231 [M] • For the generalization to functions with values in a 
Clifford algebra, called again slice hyperholomorphic or slice monogenic functions, we 
refer the reader to [HI UHl |15], and to [57] for functions with values in a real alterna- 
tive algebra. Finally, we mention that it exists a non constant coefficients differential 
operator whose kernel contains slice hyperholomorphic functions defined on suitable 
domains, [35]. The theory of slice regular functions allows to define the quaternionic 
functional calculus and its associated S'-resolvent operator. The importance of the 
S'-resolvent operator is in the definition of the quaternionic version of the operator 
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(/ — zA)~^ that appears in the reahzation function s{z) = D + zC{I — zA)~^B. It 
turns out that when A is a quaternionic matrix and p is a quaternion then (/ — pA)~^ 
has to be replaced by {I —pA){\p\'^A'^ — 2Re{p)A + I)~^ which is equal to p~^S'^^(p~\ A) 
where S]^^{p~^, A) is the right ^-resolvent operator associated to the quaternionic ma- 
trix A. For some results on the quaternionic functional calculus we refer the reader 
to [311 EH EH IH] • Slice monogenic functions admit a functional calculus for n-tuples 
of operators and for this theory we mention [l3l HHl [36] • The book [17] collects some 
of the main results on the theory of slice hyperholomorphic functions and the related 
functional calculi. 

Finally we mention the paper [T5 | [TB | [T^. where Schur multipliers were introduced and 
studied in the quaternionic setting using the Cauchy-Kovalesvkaya product and series 
of Fueter polynomials, and the papers [591 Ell |6^, which treat various aspects of a 
theory of linear systems in the quaternionic setting. Our approach is quite different 
from the methods used there. 

The paper consists of eight sections besides the introduction, and its outline is as 
follows: in Section [2] we review the main aspects of the theory of slice hyperholomorphic 
functions and the S'-resolvent operators. In Section [3] we study the counterpart of 
state space equations in the slice hyperholomorphic setting. This leads us naturally 
to the notion of rational function, defined and studied in Section jH In Section O 
we study certain positive definite functions. This paves the way to the definition of 
Schur multipliers and to a version of Schwarz' lemma in Section [61 In Section [71 we 
exhibit a coisometric realization of a Schur multiplier. Section [Hlis devoted to the Schur 
algorithm in the present setting. In the last section we discuss briefly future directions 
of research. 

2. Slice hyperholomorphic functions and the S'-resolvent operators 

In this section we introduce the preliminary results related to the theory of slice hyper- 
holomorphic functions and the quaternionic S'-resolvent operators. In the next section 
we will use these tools to define in a suitable way the quaternionic slice hyperholomor- 
phic transfer function associated to quaternionic linear systems. 

By HI we denote the algebra of real quaternions p = Xq + ixi + jx2 + kx^ which can 
also we written as p = Re(p) -|- Im(p) where Xq = Re(p) and ixi + jx2 + kx^, = lm{p). 
By S we indicate the set of unit purely imaginary quaternions, i.e. 

El = {p = Xii + X2j + x^k : xl + xl + xl = 1}. 

In the literature there are several notions of quaternion valued hyperholomorphic func- 
tions. In this paper we will consider a notion which includes power series in the 
quaternionic variable, the so-called slice regular or slice hyperholomorphic functions. 
The main reference for the material in this section is the book [17] . 

Definition 2.1. Let Q C M. be an open set and let f : Q ^ M be a real differentiable 
function. Let / G § and let fi be the restriction of f to the complex plane C/ := R + JM 
passing through 1 and I and denote by x + ly an element on C/. 
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(1) We say that f is a left slice regular function ( or slice regular or slice hyperholo- 
morphic) if, for every I G S, we have: 



d_ 

dx 



dy 



fj{x + ly) = 0. 



(2) We say that f is right slice regular function (or right slice hyperholomorphic) 
if, for every I ^ E>, we have 



—f,{x + Iy) + —fj{x + Iy)I] =0 



Definition 2.2. Given p G H, ifq is not real we can define Ip := Im(p)/|Im(j))|, so p = 
Re(p) + Ip\lm{p) I . We denote by [p] the set of all elements of the form Re{p) + J|Im(p) | 
when J varies in S. We say that [p] is the 2-sphere defined by q. 

Definition 2.3. Let Q be a domain in H. We say that Q is a slice domain (s- domain 
for short) if Q CiW is non empty and if Q (1 Cj is a domain in Cj for all I E E>. We 
say that Q is axially symmetric if, for all p G Q, the 2-sphere [p] is contained in Q. 

In the sequel we will work mainly on the unit sphere in EI with center at the origin, 
which is trivially an axially symmetric s-domain. 

Theorem 2.4 (Representation Formula). Let Q (1 M. be an axially symmetric s- 
domain. Let f be a left slice regular function on f2 C EI . Then the following equality 
holds for all p = X + Ipy G Q: 

1 

2 1' 

where z := x + ly, z := x — G H. Let f be a right slice regular function on f2 C EI. 
Then the following equality holds for all p = x + Ipy G Q: 



(2.1) 



f{p) = f{x + Ipy) = - f{z) + f{-z) 



-LI 

2 ^ 



m - m 



(2.2) 



f(x + Ipy) = -[fiz) + fiz) 



The Representation Formula allows to extend any function / : f2 C C/ — t- H defined 
on an axially symmetric open set Q intersecting the real axis and in the kernel of the 
Cauchy-Riemann operator to a function / : f2 C H — )■ EI slice regular where Q is 
the smallest axially symmetric open set in EI containing Q by means of the extension 
operator 



(2.3) ext(/)(p):=-[/(z) + /(z) 



1 



z, z E Cj. 



Slice regular functions satisfy an identity principle, specifically, two of them coincide 
in a domain Q intersecting the real axis if they coincide on a subset of fl M having 
an accumulation point. As a consequence we have: 

Proposition 2.5. Any quaternion valued real analytic function defined in (a, 6) C M 
extends uniquely to a slice regular function defined in a suitable open set containing 
{a,b). 

Proof. Let / : (a, 6) C M — )■ EI be a real analytic function. Then it can be uniquely 
extended to a suitable domain D containing {a,b), by extending its real components, 
to a function / : D C C ^ M which is in the kernel of the Cauchy-Riemann operator. 
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By using the extension operator fl2.3p / extends to a slice hyperholomorphic function 
defined on the smallest axially symmetric open set Q containing D. The uniqueness 
follows from the identity principle. □ 

Given two slice regular functions their product is not, in general, slice regular. It 
is possible to introduce a suitable product denoted by ic, see [171 Definition 4.3.5, p. 
125] (note that in the hterature this product is denoted by * but here we use this 
symbol in order to avoid confusion with adjoint of operators). Here we describe the 
T*r- product which gives a left slice regular product. It is possible to define an analog 
product to multiply two right slice regular functions and obtain a function with the 
same regularity. Let f2 C EI be an axially symmetric s-domain and let f,g : f2 — EI 
be slice regular functions. Fix /, J G S, with I A. J. The Splitting Lemma guarantees 
the existence of four holomorphic functions F,G, H, K : f2 fl C/ — C/ such that for 
all z = X + ly e n Cj 

fj{z) = F{z) + G{z)J, gj{z) = H{z) + K{z)J. 

Define the function fi gj : f2 fl C/ — ?■ H as 

(2.4) /, ^ gj{z) = [F{z)H{z) - G{z)K{i)] + [F{z)K{z) + G{z)lm]J. 

Then fj-k gj[z) is obviously a holomorphic map and hence we can give the following 
definition. 

Definition 2.6. Let Q C M be an axially symmetric s-domain and let f,g : f2 — )> EI 
be slice regular. The function 

f-kg{p) = ext{fj-k gi){p) 
is called the slice regular product of f and g. 

Remark 2.7. It is immediate to verify that the ^-product is associative, distributive 
but, in general, it is not commutative. 

When slice regular functions are defined on a ball C EI with center at the origin 
and radius R then they admit power series expansions. Suppose that on Br we have 
f{p) = J2n=oP''^n and g{p) = Xl^oP"^™, then the ^-product is given by 

oo n 

(2.5) (/*fl')(p) = X]^"^"' c„ = ^a^6„_^. 

n=0 r=0 

In other words, the coefficient sequence {Cn)nmo the Cauchy product, or the convo- 
lution product, of the sequences of coefficients associated to / and g. See for instance 
(2) p. 199] for the Cauchy product in a non -commutative setting. 

Remark 2.8. From Proposition 12. 5[ it immediately follows that f -k g is uniquely 
determined by f{x)g{x), x G Mr fl M. 

We observe that another product in Clifford analysis (see [271112]), namely the Cauchy- 
Kovaleskaya product which allows to obtain Cauchy-Fueter regular functions or mono- 
genic functions, can be seen as a (different) convolution. More generally, pointwise 
product is often best replaced by convolution of underlying coefficient sequences. See 
[27] . We mention as an example the Wick product in white noise space analysis. See 
[60] , and see [lOl [11] for applications to linear system theory. 
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Remark 2.9. Let f{p) = X^^p^a^ and g(j)) = J2nP"'^n- If f has real coefficients then 

f-*^9 = g-^^f- 

Given a slice regular function / it is possible to construct its slice regular reciprocal, 
denoted by /~*. We do not provide the general construction, which can be found in 
[17] since it is sufficient to construct the inverse of a polynomial or a power series with 
center at the origin. Given /(g) as above, let us introduce the notation 

oo oo n 

rip) = XI P"^"' f'iP) = = ^P^Cn, Cn = ^arttn-r- 

n=0 n=0 r=0 

Note that the series /** has real coefficients. The slice regular reciprocal is then defined 
as 

In an analogous way on can define the right slice regular reciprocal /~* := f'^{f^)~^, 
of a right regular function f{p) = Ylm'^nP^- 

We do not introduce a specific symbol in order to distinguish the left and the right slice 
regular reciprocal, since the context will clarify the case we are considering. There is 
however a remarkable case of reciprocal that deserves further explanations. Consider 
the function 

S{r, p) := p — q, p,q & EI. 
Its slice regular reciprocal can be constructed in four possible ways: we can construct 
a reciprocal which is left (resp. right) regular with respect to the variable p or left 
(resp. right) regular with respect to the variable q. Accordingly to these possibihties 
we obtain the function (see [17]) 

Sl\p, q) = -iq' - 2qRe{p) + \p\Y\q - p) 

which corresponds to the left slice regular reciprocal in the variable q. The function 
S]^^{p,q) is left regular in the variable q by construction and it turns out to be right 
regular with respect to p. While the other possibility gives 

Sn\p,q) := -{q - p){q' - 2Re{p)q + \p\^)-\ 

which is right slice regular in q and left slice regular with respect to p. 

Remark 2.10. When no confusion will arise we will write instead of {p—q)~* its explicit 
expression using S~[^{p,q) or S]^^{p,q), according to the left or right slice regularity 
required. 

It is possible to show that both kernels can be written in two different ways 

Sl\p, q) = -iq' - 2gRe(p) + \p\Y\q -P) = ip- q)ir' - 2Re(g)p + \q\Y\ 

and 

Sr\p, q) = -{q - P){q^ - 2Re(p)g + \p\^Y^ = {p" - 2Re(g)r + \q\^r\r - g), 
thus we have 

SR\p,q) = -S2\q,p). 
Given two polynomials P{p), Qip) it is possible to define two possible "quotients" as 
P~* -kQ oi Q -k P~* and due to the noncommutativity of quaternions, these quotients 
do not coincide. Thus it is necessary to make a choice between the left and the right 
quotient of P and Q and this will be one of the various definitions of rational functions. 
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A number of equivalent characterization of rational functions are given in Section HI 

The realization of Schur functions in the form s{z) = D + zC{I — zA)~^B, where one 
wishes to replace the complex variable 2; by a quaternionic variable, and where the 
operator matrix 

'A 

is now a quaternionic operator matrix, requires a new concept to replace the classical 
resolvent operator (/ — zA)~^ . This new object is the so called S'-resolvent operator 
associated to the quaternionic functional calculus. Let ^ be a two sided quaternionic 
Banach space, we denote by BiV) the quaternionic Banach spaces of all bounded linear 
(left or right) operators endowed with the natural norm. A remarkable fact is that 
the following theorem holds in the case in which the components of the quaternionic 
operator A do not commute. 

Theorem 2.11. Let A e BiV). Then, for \\A\\ < \r\, we have 
00 

(2.6) J2 ^V-i-" = -{A^ - 2Re{r)A + \r\^I)-\A - rl), 

n=0 

(2.7) ^r-^-M" = -{A - Ir){A^ - 2Re{r)A + \r\^I)-\ 

n=0 

where I denotes the identity operator. 

The notion of S'-spectrum of a quaternionic operator A is suggested by the definition 
of S'-resolvent operator that is the kernel for the quaternionic functional calculus; see 
dM} and fET|) . 

Definition 2.12 (The S'-spectrum and the S'-resolvent sets of quaternionic operators). 
Let T G B{V). We define the S-spectrum as{A) ofT as: 

as{A) = {r G H : - 2 Re{r)A + |r|^/ is not invertible}. 

The S-resolvent set ps{A) is defined by 

ps{A)=m\as{A). 

Remark 2.13. We have proved that the S'-spectrum as{A) is a compact nonempty 
set and if p G EI belongs to as (A), then all the elements of the sphere [p] are contained 
in as{A). 

Definition 2.14 (The S'-resolvent operator). Let V be a bilateral quaternionic Banach 
space, A G B{V) and r G psi^)- We define the left S-resolvent operator as 

(2.8) S2\r, A) := -{A^ - 2Re{r)A + \r\^iy\A - rl). 
We define the right S-resolvent operator as 

(2.9) S^\r, A) := -{A - If){A^ - 2Re{r)A + |rp/)-i. 



SCHUR FUNCTIONS IN THE SLICE HYPERHOLOMORPHIC SETTING 9 

Theorem 2.15 (The S'-resolvent equation). Let V be a bilateral quaternionic Banach 
space, A G BiV) and r G ps{A). Then the left S -resolvent operator defined in h2. <^)) 
satisfies the equation 

(2.10) Sl\r,A)r-ASl\r,A)=I. 

and the right S-resolvent operator defined in 1^2. 9\) satisfies the equation 

(2.11) sS^\s,A) - S^\s,A)A = I. 

Moreover, S2^{r, ■) : ps{A) — )■ B{V) is right slice regular in the variable r and S]^^{r, ■) : 
PsiA) — )■ B(y) is left slice regular in the variable r. 

The notation introduced in Remark 12.101 is not necessarily valid when we replacing 
operators in place of the quaternionic variables. We prove that it is still the case in 
the following proposition which is crucial and can be proved by suitably modifying the 
proof of Theorem 12.111 

Proposition 2.16. Let be a two sided quaternionic Hilbert space and let A be a 
bounded right linear quaternionic operator from H into itself. Then, for \p\ \\A\\ < 1 
we have 

oo 

(2.12) ^pM" =p-i5^i(p-\A) = {I-pA){\p\'A'-2Re{p)A + ir\ 

n=0 

and 

oo 

(2.13) (J-M)-* = ^pM". 

n=0 

Proof. From (12. 7p . it is immediate that 

oo 

Y,p"A"=p~'S^\p-\A) 

n=0 

and by writing explicitly S]^^{p~^,A) we obtain 

and with some computations we get the second equality in fl2.12p . To prove (I2.13P we 
consider the function p~^ — q and its slice regular reciprocal with respect to q which 
is obtained using the formula /~* = /'^(/*)~^, see Section 2. We obtain {p~^ — *?)"* = 
S]^^{p~^, q) and, by the functional calculus, we can substitute g by a quaternionic oper- 
ator A (it is crucial to observe that the components of A do not necessarily commute) 
and we obtain: 

{p~H - AY' = -S^^{p~\A) = -{A-p~H){A^ - 2Re(p-i)A + b"^/)-^ 

since 

p-\p-^I - Ay = [p{p-^I - A)]-'' = (/ - pA)-\ 
By using f l232|) we obtain (|27[3|) . □ 
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3. State space equations and realization 
We now show that if we consider the quaternionic hnear system 

( Xn+l = XnA + UnB, n = 0, 1,... 

where A, B, C, D are given quaternion matrices, then the quaternionic "transfer 
function" cannot be defined by simply taking 2^-transform as in the complex case. The 
following considerations show the problem that arises: Given a sequence of quaternions 
(or of quaternionic matrices of common size) U = {un}neN, we define the quaternionic 
-E-transform as 

oo 

Z{U) ■.= U{p) :=5^p"u„, 

n=0 

SO the 2^-transform is right linear 

Z{UA) = Z{U)A. 

Since Z{U) is a power series centered at the origin it is slice hyperholomorphic. The 
main properties of the classical ^([/)-transform still hold in the quaternionic setting. 
In particular, if we set 

r_if/ := (mi,M2, . . .) 

we have 

Z{t^iU) = p~^Z{U) if Mo = 0. 
If we apply the ^-transform to system (13. ip we get 

p~^X{p) = a:{p)a + u{p)b, n = 0,l,... 

y{p) = x{p)c + u{p)D. 

The natural definition of the transfer function of the system is 

n{p) := mp)r'y{p) 

so we have to solve the quaternionic equation 

p^^X{p) - X{p)A = BU{p) 

which has the solution 

oo 

X{p) = Y^p''BU{p)A''+\ 

n=0 

Because of the term B, this expression need not be slice hyperholomorphic. By replac- 
ing this solution in the second equation in (13. 2p . we obtain 

oo 

y{p) = J2p''BU{p)A''+^C + U{p)D. 

n=0 

It turns out that 'H(p) depends on U{p), in fact 

oo 

n{p) := {U{p)Y\'Y,P""BU{p)A''+^C + U{p)D 



(3.2) 



n=0 

oo 



(W(p))"^ ^p"5W(p)A"+iC + D, 



n=0 
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and need not be slice hyperholomorphic. We now follow an approach based on slice 
hyperholomorphic to overcome the above difficulty and to give a good definition of 
transfer function. 

Theorem 3.1. Let A, B, C, D be given matrices of appropriate size with quaternionic 
entries. Suppose that {wnjnGNo given sequence of vectors with quaternionic entries, 
and of appropriate size. Consider the system 



(3.3) 



Xn+l = AXn + BUn, n = 0,l, 
Un = CXn + DUn, 



and define its slice hyperholomorphic transfer function matrix-valued as 

n{p) := yip) ^ (w(p))-* 

where y{p) andU{p) are the slice hyperholomorphic extensions of the Z-transforms of 
yn and of Un, respectively. Then we have 

(3.4) n{p) = D+pCi.{I -pAy*B. 

Moreover, in terms of the right S-resolvent operator Iji3.4\ ) can be written as 

(3.5) H{p) = C i.Sj^Xp-\A)B + D. 

Proof. Let us first consider system (13. 3 P in the complex plane C/, for a fixed / G S, 
where A, B, C, D, now denoted by a, b, c, d, are given element in C/. Suppose that 
{Un}neNo is a given sequence of complex numbers in Cj: 

, . ( Xn+l = aXn + bUn, 72 = 0,1,... 

\ yn = CXn + dUn- 

By setting 

z := u + Iv E Cj 
and by taking the classical 2^-transform we get 

/ X{z)=za X{z) + zbU{z) 
^'^■'^ \y{z) = cX{z) + dU{z). 

On C/ all the objects commute, thus we have: 

X{z) = (l-za)-^ zbU{z), 



(3. 



y{z) = c X{z) + dU{z). 



All the functions in the complex variable z involved in system (13. 8p are holomorphic 
on the plane C/, thus we can use the extension operator (12. 3p to obtain a function in 
the quaternionic variable p which is slice hyperholomorphic and since 



we get 
(3.9) 

from which 
(3.10) 



e^i{f{z)g{z)) = ext{f{z)-kg{z)) = f{p)-kg{p) 



M[p) = {l-pay*i<{pb)i<U{p), 
y{p) = c * X{p) + d U{p) 



X{p) = {l-pa)-*^{pb)^U{p), 

y[p) =c i<{l - pay* -k{pb) -kU{p) + d -k U{p). 



12 
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Since we have defined 'H(p) as 
we obtain 

n{p) = (c * {I -pay* * (pb) i.u{p) + d u{p)^ * {u{v)y* 

= c * (1 — pay* (pb) + d. 

This function is slice hyperholomorphic in p with values in the quaternions so it makes 
sense to consider a, 6, c and d as quaternions. In order to get a matrix valued function 
it is sufficient to replace a, 6, c and c?, respectively, with the matrices A, B, C, D, of 
appropriate size and with quaternionic entries. We obtain the quaternionic transfer 
function 

n{p) = D + c ^ (/ - pAy* * {pB) = D + pc * (/ - pAy*B 

which is a slice regular function. Finally, using (12.121) . we obtain the equality (13. 5p . □ 
Remark 3.2. Using the explicit expression for (1 — pA)~* we can also write 

H{p) = D + C ^ (J - pA){\p\^A^ - 2Re{p)A + J)"^ ^ (pB), 
and if we remove the ^-product we get 

nip) = D + {pC- \p\^CA){\p\''A'' - 2Re{p)A + ly^B. 

4. Rational functions 

Motivated by the discussion in the previous section, we now give various equivalent 
definitions of rationality. We first consider the case of H*^^ ^-valued functions of a real 
variable, and give a definition and some preliminary lemmas. 

Definition 4.1. A function f{x) of the real variable x and with values in H*-'^^^ will 
be said to be rational if it can be written as 

(4.1) /(.) - 



m[X) 



where M is a M -valued polynomial and m E M[x]. 

Clearly sums and products of rational functions of appropriate sizes stay rational. The 
case of the inverse is considered in the next lemma. 

Lemma 4.2. Let f be a rational function from M into H^^^, and assume that /(xq) 
is invertible for some Xq G M. Then, f{x) is invertible for all x E at the possible 
exception of a finite number of values, and f^^ is a rational function. 

Proof. We proceed by induction. When = 1, a polynomial in M.[x] can be seen as a 
matrix polynomial in M^^^[x], and this latter is invertible for all but at most a finite 
number of real values of x. Let now r{x) = ^j^, where ti G M.[x] and t G ]R[x]. Assume 
t ^ 0. Then, since x is real, we may write 



tix) t{x) 



t{x)t{x) {tt){xy 
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This ends the proof for = 1 since {tt){x) G M.[x]. 

Assume now the induction proved at rank and let i? be a H'-^"'"^^^^^"^"'^-* -valued 
rational function, invertible for at least one xq G M. Then identifying M*^^"*"^''^*^^"'"^'' 
with ]R4(Af+i)x4(Af+i) that R is invertible for all, but at most a finite number of 

real values of x 

^a(x) b{x) 
^c(x) d{x) 

Without loss of generality, we can assume that a{x) is not identically equal to (oth- 
erwise, multiply R on the left or on the right by a permutation matrix; this does not 
change the property of R or of R~^ being rational). We write (see for instance [521 
(0.3), p. 3]) 

a{x) b{x)\ _ f I 
c{x) d{x) J \c{x)a{x)~^ In 



R{x) 



X 



X 

a{x) \ fl a{x)-%x] 

d{x) - c{x)a{x)-^b{x) ) \0 /„ 



and so d{x) — c{x)a{x) ^b{x) is invertible for all, but at most a finite number of, values 
X G M. We have 



R-Hx) 



(4.2) 



-a{x) ^b{x) 



X 



a(x)-^ 



{d{x) — c{x)a{x) ^b{x)) \ —c{x)a{x) 







The induction hypothesis at rank implies that {d{x) —c{x)a{x) ^b{x)) ^ is rational, 
and so is R~^, as seen from (14. 2p . □ 

Lemma 4.3. Consider a polynomial M{p) G H[^^^[p].- 



Then 

where D = Mq, 



m{p) = J2p^Mj. 

j=0 

M{p) =D+pC*{I- pA)-*B 

\ 



/Oat In ■ ■ ■ 

On On In On 



B 



A 



/On\ 

On 

\InJ 



On ■■■ 

\0n On 



On In 
On On/ 



C = {Mj Mj.i ■■■ Ml) . 
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Proof. The equality easily follows from the fact that 

/ In pIn P^In ■ ■ ■ P'^In \ 

On In pIn ■■■ P'^^^In 

{I-pA)-^= : : 

Otv In pIn 

\0n On ■■■ On In J 



□ 



Assume now that the H^^^-valued function f{p) can be written as 

(4.3) f{p)=D+pC^{I-pAy*B, 

where A, 5, C and D are matrices with entries in H and of appropriate sizes, and that 
D is invert ible. Then, 

(4.4) f{pY^ = D-^ - pD-^Ci. (/ - p{A - BD-^C))-*B 

We also recall the following formula. See for instance [23] for the case of rational 
functions of a complex variable. The proof is as in the classical case, and will be 
omitted. 

Lemma 4.4. Let 

f,{p) = D, + pC, ^ {In, - pA,)-*B,, j = 1, 2 

he two functions admitting realizations of the form fl4.3p . and respectively H*^^^ and 
M^^^-valued. Then the M'^^^^ -valued function /i^/2 can be written in the form (14.31) . 
with D = D1D2 and 

The proof of the following Lemma is immediate. 
Lemma 4.5. Let 

f,{p) = D, + pC, ^ {In, - pA,)-*B^, j = 1, 2 

be two functions admitting realizations of the form (14. 3p . and respectively H*^^^ and 
M^^^-valued. Then the M.^'^^^^^^'^ -valued function (/i f2) can be written in the form 
K^ . with D = [Di D2) and 

In the following theorem, 

00 

/(p) = ^diag(p",p",...,p")/„, 

M times 

which we will denote for short by 



n=0 
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is a series of the quaternionic variable p, with coefficients in H*^^^, converging in a 
neighborhood of the origin. We say that / is rational if any of the five conditions listed 
in the theorem holds. 

Theorem 4.6. Let f be a M'^^^ -valued function, slice hyperholomorphic in a neigh- 
borhood of the origin. Then, the following are equivalent: 

(1) There exist matrices A,B and C, of appropriate dimensions, and such that 

(4.5) fn = CA^-^B, n = l,2,... fo = D. 

(2) / can be written as (14. 3p 

f{p) =D+pCi.{I- pA)-*B. 

(3) The right linear span Ai{f) of the columns of the functions Rof, Rof, . . . is a finite 
dimensional right quaternionic Hilbert space. 

(4) The function f{x) is a rational function from M into H*''^^^. 

(5) The entries of f are of the form P -k Q~* , where P and Q are slice holomorphic 
polynomials, and Q{0) ^ 0. 

Proof. In the proof we adapt well known arguments from the theory of matrix-valued 
rational functions to the present case. For similar arguments in the classical case, see 
for instance |24j . 

Assume that (1) is in force, and set D = sq- Then, for p such that \p\ ■ \\A\\ < 1, the 
series ^^o^*"^" converges in M^^^, and 

oo 

= (/ - pA)-\ where I = In, 

n=0 

and we can write 

oo 

J2p''fn = D+pC*iI-pA)-*B, 

n=0 

SO that (2) is in force. 

Assume now (2). Then, 

R^f = Ci.{I -pAy^A'^-^B, n = l,2,..., 

so that Ai{f) is spanned by the columns of D and the columns of the function 
C -k {I — pA)~*, and is in particular finite dimensional, so that (3) is in force. 

Assume now (3). Then there exists an integer mo G N such that for every m G N and 
V there exist vectors mi, . . . , Umo such that 

(4.6) R^^fv = Y,R^fum. 

m=l 

Of course, the Uj need not be unique. Let E denote the H^^™" "^-valued slice hyperholo- 
morphic function 

E={Rof Rlf ... i?r/)- 
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Then, in view of (14. 6p . there exists a matrix A G H'"°''^™°'^ such that 

RqE = EA, 



so that 
and so 
and 



E{p) - E{0) = p ^ E{p)A = E{p) * pA, 
E{p) = E{0) {I - pA)-* 



{R,f){p) = E{p) 







Thus we have 



f{p)-f{Q)=pE{Q)^{I-pA)- 







Vo/ 



It follows that / is of the form f l4.3p . and so (2) holds. Since (2) implies trivially (1), we 
have shown that the first three claims are equivalent. We now turn to the equivalence 
with the other two characterizations. 

Assume that (3) holds. Then, the restriction of f l4.3p to p = x G M gives 

f[x) = D + xC{I -xAy^B. 

and thus / is a rational function of x in view of Lemma 14.21 that is (4) is in force. 
Assume now (4), and let f{x) = as in (14. ip . Since / is assumed defined the origin, 
we may assume that m(0) ^ 0, and f{x) is the restriction to the real line of the slice 
hyperholomorphic function 

f{p) = M{p) -k m{p)~* = m{p)~* -k M{p). 

Note that since m has real coefficients we also have m{p)~* -k M[p) = m{p)~^M{p). So 
(5) holds. Assume now (5) in force. Then, in view of Lemma 14. 4[ each of the entries 
of / admits a realization, and the result follows since the property of realization stays 
under concatenation. Use Lemma 14.51 and its counterpart for columns rather than 
rows. □ 

We note that a short introduction to realization theory in forms of a sequence of 
exercises can be found in [3l pp. 328-329]. 

5. Positive definite functions and reproducing kernels 

The H-valued function k{u,v) defined for u,v in some set Q (in the present paper, Q 
will be most of the time the unit ball of H) is called positive definite if it is Hermitian: 



(5.1) 



k{u,v) = k{v,u), Wu.v^Q, 
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and if for every N E every Ui, . . . , Uat G and Ci, . . . , Ctv G HI it holds that 

N 

C£k{u£, Uj)cj > 0. 

(Note that the above sum is a real number in view of (15. ip ). As in the complex- valued 
case, associated to /c is a uniquely defined reproducing kernel quaternionic (right)- 
Hilbert space Jif{k), with reproducing kernel k{u,v), meaning that: 

(1) The function u i— )• k{u,v)c belongs to Jf{k) for every choice of v E fl and c G H, 
and 

(2) it holds that 

for every choice of / G Jif{k) and of m G and c G H. See for instance Theorem 
8.4, p. 456] for a proof and [T3l Proposition 9.4, p. 458] for a proof of a result that, 
with similar techniques, proves also the following lemma. 

Lemma 5.1. Let Jif{k) be a right reproducing kernel quaternionic Hilbert space ofM- 
valued functions defined on a set fl, with reproducing kernel k{u, v) . Then, the function 
f belongs to J^{k) if and only if there exists M > such that the kernel 



The smallest such M is equal to the norm of f in M'{k). 

Reproducing kernels are a main tool in operator theory. Reproducing kernels of the 
form c{zw), z,w E where c{t) = ^^^c„t", c„ > 0, for all n G N is an analytic 
function is a neighborhood of the origin are an important case. In the case of Hardy 
space, Bergman space, and Dirichlet space the functions c(t) are given by 

cit) = ^ c(t) = c(t) = - ln(l - 1). 

We will consider the generalization of these kernels to the case in which the variables 
considered are quaternions, i.e. we will consider 

+ 00 

(5.2) c(M) = 5Zc„p"g", q,pem. 

n=0 

Note that due to the noncommutative nature of quaternions, the series in (15. 2 p does 
not coincide with ^n=o '^ri{pq)^ ■ 

Theorem 5.2. Let A G Nq and let {cn)n£A be a sequence of strictly positive numbers, 
and assume that 

R = i ^ > 0. 

hm sup^g^ Cn 

Then, the function 
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is positive definite in the ball \p\ < R. The associated reproducing kernel Hilbert space 
consists of the functions 

f{p) = Y.p^f-^ 

with coefficients /„ G H such that 



< oo. 



In the case of the Hardy space we have 

Proposition 5.3. The sum of the series Yln=oP^^ ^'^ function k{p,q) given by 
(5.3) k{p,q) = {1 - 2Re{q)p + \q\yy\l - pq) = {1 ~ pq){l - 2Re{p)q + \p\^f)-\ 

The kernel k{p,q) is defined for p ^ [q~^] or, equivalently, for q ^ [p''^]- Moreover: 

a) k{j), q) is left slice regular in p and right slice regular in q; 

b) k{p,q) = k{q,p). 

Proof. Consider the function 

kq{z) = z = x + Iy. 
1 — zq 

The proof of the first equality is an application of the representation formula (12. ip . 
Consider now the function 

kp{w) = 

1 ~ pw 

the representation formula (12. 2p gives k{p,q) = (1 — pq){l — 2Re(p)g+ \p\'^q'^)"^- The 
function (1 — pq){l — 2Re(p)g + \p\'^q'^)~^ is right slice regular in the variable q by 
construction and it is slice regular in the variable p in its domain of definition, since it 
is the product of a slice regular function and a polynomial with real coefficients. By 
the identity principle it coincides with the first expression which is slice regular in p by 
construction. Assertion a) follows. Point b) follows from the chain of equalities 



k{q,p) = [{1 - pq){l - 2Re{p)q + \p\^^)-^ 
= (1 - 2Re(p)g + \p\'^f)~\l - qp) = k{q,p). 

□ 

We note (see for instance [501 12]): 
Theorem 5.4. The following are equivalent: 

(1) The function s is analytic and contractive in the open unit disk. 

(2) The function s is defined in D and the operator of multiplication by s is a contraction 
from the Hardy space i?2(D) into itself. 

(3) The function s is defined in D and the kernel 



1 - s{z)s{w) 

ks{z,w) = — = — 

1 — zw 

is positive in the open unit disk. 
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We prove in Theorem 16.21 the counterpart of the above theorem in the quaternionic 
setting. 

Note that kg can be rewritten as 

oo 

This is the form which will be used in the quaternionic setting. See (16. 3 p in Section [6l 

The reproducing kernel Hilbert space with reproducing kernel ks{z,w) was first intro- 
duced and studied by de Branges and Rovnyak; see [291 EH] . This space will be denoted 
by Jif{s). It is equal to the operator range ran \/I — MgM*, where Ms denotes the 
operator of multiplication by s from H2{Vi) into itself, that is, Mg is the commutative 
version of the operator defined by (16. 2p in the following section. 

The space Jif{s) is the state space for a coisometric realization of s. Furthermore, s is 
a Schur function of and only if the function defined by 



(5.4) s\z) = s{z) 

is a Schur function. The space J^{s'^) is the state space for an isometric realization of 
s. A unitary realization for s is given in terms of the reproducing kernel Hilbert space 
with reproducing kernel 



(5.5) Ds{z,w) 



z — w 

\ = kst{z,w) I 

\ z — w / 



For more information on the Schur algorithm, see for instance the papers [61] and the 
books [2l[l9lll8]. 

6. Schur multipliers and Schwarz' lemma 

We will define a finite energy quaternionic signal as a sequence x = {xo,Xi,...) of 
quaternions such that 



oo 



I II 2 dof . V ^ I 1 2 

\X\\ = > \Xr,\ < OO 



n=0 

that is, an element of £2(1^0) H)- Note that £2(^0,13) is a quaternionic right Hilbert 
space when endowed with the H-valued inner product 

00 

[2;, y] = ^y^Xn, where y = {yo, yi, . . .). 

n=0 

To X G ^2 (No, H) we associate the function 

00 

fx{p) = XI P"^'^ 
n=0 
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of the quaternionic variable p. We denote by H2 the right quaternionic vector space of 
such power series, endowed with the inner product 

00 

[fx: fy] — ^ ^ UnXn- 
n=0 

To ease the notation, we will use the notation x{p) rather than fx{p)'- 

00 

(6.1) x(p) = ^p"x„. 

n=0 

Definition 6.1. We set £(No,lHI) denote the set of sequences of quaternions, indexed 
by No, and io{^o, EI) denote the set of sequences of i{No, M) for which at most a finite 
number of elements are non-zeros. 

The shce regular product of x G £2(^0, H) and y G £9(^0, H) is defined by 

n 

(x -k y^n ^ ^ ■^myn—mi 0,1,... 
m=0 

Note that 

xi^l = x, where 1 = (1, 0, 0, . . .). 
We now state the counterpart of Theorem 15.41 in the quaternionic setting: 

Theorem 6.2. Let s G £(No,H) (see Definition \6. Then, the following are equiva- 
lent: 

(1) The map x ^ si^ x, first defined for x G ^o(I^O)EI), extends to a contraction from 
4(No,e) into itself. 

(2) The series s{p) = Yl'^=QP"'Sn converges in M and the operator 

00 

(6.2) (M,x)(p) = ^p"s(p)x„ 

n=0 

is a contraction from the quaternionic Hardy space H9(B) into itself, where x{p) = 

Er=oP"^neH2(B). 

(3) The series s{p) = Xl^oP"-^" converges in M, and the function 

00 

(6.3) h{p, q) = ^p"(l - s(p)^)g" = (1 - s{p)J{^) * (1 - pq)"* 

n=0 

is positive definite in the unit ball B of H. 

Proof. We first assume that (1) is in force. Since s -k 1 = s we note that in fact 
s G £2(^0, H). The power series s{p) makes thus sense for every p in the open unit 
ball of H. We denote by Ms the image of the operator x ^-^ s-kx under the map (16. ip . 
Then, for x G i'o(No, H) we have 

00 00 00 00 

(6.4) J]p"5(p)a:„ = ^^^"(J^ p"^3„.)a:„ = ^^^"(s^a;)^, 

n=0 n=0 m=0 n=0 

where the sum is finite since x G £o(f^05EI). Since in a reproducing kernel Hilbert 
space, convergence in norm implies pointwise convergence, we see that the operator 
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Ms is given by (16 ■4p for every x G £2(1^0) U)? ^i-iid hence (2) is in force. 

Assume now that (2) holds. Consider the kernel k{p,q) in (15.31) . which on \pq\ < 1 
admits the power series expansion 

n=0 

For sake of simplicity we will write kq{p) instead of k{p, q). Let qi,q2 and Ci, C2 G EI 
and compute 

c^{M*kg^ci) (^2) = [M*kg,ci , /cg2C2]H2(B) 
= [kg^ci , Ms(/c52C2)]h2(b) 



[kq.ci, Ep"s(p)g2"c2]H2 



n=0 



[^pns(p)q^C2, kq,Cl]u2 
n=0 

00 

C2 1 Yl ^2Mr"ci 



vn=0 

Thus, the adjoint of Mg is given by the formula: 



(M;(fc(g))(p) = X^p".(g)"r 



n=0 

and (3) is just a rewriting that the operator / — MgM* is non negative. 

Assume now that (3) holds. We consider the relation (that is, the linear subspace) 
in H2(B) X H2(B) spanned by the pairs 

J2 P^<fc , p''s{p)4q)^c 

,n=0 n.=0 

where p runs in IB and c in H. The domain of Rg is dense, and the positivity of the 
kernel implies that Rg is a contraction, meaning that if {f,g) G Rs, then 



||fi'||H2(B) < Hi llHai , 

It follows that Rs extends to the graph of an everywhere defined contraction, which we 
will denote by T. We now compute T*: 

C^{T*kq^Ci){qi) = [T*kq^Ci, kq^C2\u2(M) 

00 

= [kq.Ci, Ep"s(g2)g2''C2]H2(B) 
n=0 

(00 
X]^2s(g2)grci 
n=0 

= c^{Msikq,c^))iqi). 
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Thus T* = Mg. We obtain (1) by looking at the operator on each p"-. □ 
We note that the kernel fc, can also be written as 



ksij), g) = (1 - 2Re(g)p + |g| V) ^(1 - pq) * (1 - s{p)s{q)) 

= (1 - 2Re{q)p + \q\^p^)~^{l - pq - s{p)s{q) + ps{p)qs{q)). 

The operator Mg is the non-commutative version of the operator of multiplication by 
s. We note the following properties of the operators M^: 

Proposition 6.3. Let si, S2 and s be Schur multipliers depending on the quaternionic 
variable p. Then: 

(6.5) Ms.Mg, = M,,,s2, 

(6.6) MsMp = MpM, = Mps. 

Proof. The first equality follows from the associativity of the ^-product. The second 
property is a consequence of Remark 12.91 while the last one follows from the previous 
ones. □ 

The proof of the following proposition follows the case of analytic functions in the 
open unit disk (see [631 Lemma 1, p. 301] for instance for the latter), and is given 
for completeness. In the case of analytic functions in a half-plane, the corresponding 
result is more difficult to prove, and is called the Bochner-Chandrasekharan theorem. 
See [251 El]. 

Remark 6.4. We now discuss a feature which arises because we are working in a 
non-commutative setting. When considering a vector space on H, the set of right 
linear operators acting on V has a structure of vector space over EI only if \^ is a two 
sided vector space, in fact right linear operators form a left vector space, see Section 
2. However, if in V it is defined a multiplication * among vectors compatible with the 
right vector space structure, or at least, if for some w E V the multiplication w v is 
an operation in V for all v G V, then this multiplication allows to define right linear 
operators which form a right linear vector space. Thus, in this case, it is not necessary 
to require that V is two sided vector space. Indeed, we can define a right linear operator 

M^(f ) := w -kv. 

The linearity is obvious and the structure of right linear space on the operators is given 
by 

{M^X){v) ■.= wX^v, AeM. 

In particular, if = H2(]B) and M is the multiplication on the left by s{p) G V 
(assuming that the multiplication is defined for every element in V^) then we have 

Ms{f){p) = s{p)'kf{p) and if 7r{p) = ao+pai^ hp"a„ then M^(7r)(p) = s{p)'k7r{p) = 

s{p)ao + ps{p)ai + . . .+ p'^s{p)an. 

Proposition 6.5. Let T be a right-linear contraction from H2(B) into itself, which 
commutes with Mp. Then T = Ms for some Schur multiplier s. 
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Proof. Denote s{p) = {Tl){p). Then, for a G H C H2(B) (in other words, in (16. 7p and 
(16. 8p below we view a G EI as a function in H2(B) and, by the assumed right hnearity, 
we have Ta = (Tl)a) 

(6.7) ((TM,)(a)) (p) = (M,(Ta))(p) = {{Mp{Tl))a){p) = ps{p)a. 
More generally, an induction shows that 

(6.8) ((rMpn)(a))(p) = (Mpn(Ta))(p) =p"s(p)a, n = 2,3,... 
Thus, for every polynomial 

7r(p) = ao+pai^ hp"a„, 

we have 

n 

(6.9) {Tn){p) = Y,P^s{p)a,. 

j=0 

Let now / G H2(]B), with power series expansion 

oo 
3=0 

and let /„ be the polynomial 

n 

fn{p) = X^P^Oj. 
j=0 

We have lim„_j.oo ||/ — /n||H2(B) = 0, and so, by continuity of T, 

lim \\Tf - T/„||h2(b) = 

In a reproducing kernel Hilbert space convergence in norm implies pointwise conver- 
gence. Thus we have that, for every p G B, 

oo 

i.Tf){p) = ^p^s{p)aj, 

where we have used ( 16.90 . Thus, T = Mg. Finally, s is a Schur multiplier since T is 
assumed contractive. □ 

The following result is a counterpart of Schwarz's lemma for Schur multipliers. 

Theorem 6.6. Let s be a Schur multiplier, and assume that s{0) = 0. Set s{p) = 
ps^^\p). Then s^^^ is also a Schur multiplier. 

Proof. Let ,y^{s) be the reproducing kernel quaternionic Hilbert space with reproducing 
kernel ks{p,q). Since s(0) = we have that 1 = ks{p,0) G =^(s), and = 
ks{0,0) = 1. Hence the kernel 

ks{p,q) - 1 

is positive definite in B; see Lemma 15. 1[ Thus 



oo 

J2p" (l -ps(i)(p)7(%)g) r > 1, 



n=0 
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which can be rewritten as 

p ( X^P" (l - j g > 0. 

\n=0 / 

Thus, 

oo 

5^p"(i-.«(p)i(%))r >o, 

and s^^-* is a Schur multipher. □ 

We conclude this section with the following characterization of the space J^{s): 

Theorem 6.7. Let s he a Schur multiplier. Then the associated reproducing kernel 
Hilbert space J€'{s) is equal to the operator range ran 

(J-M,M;)1/2 endowed with the 

norm 

II (/ - m,m;)V2«|u.(,) = ||(j - n)u\\nm^ 

where ir is the orthogonal projection on the kernel of (/ — M^M*)^/^. 

The proof is as in the classical case and will be omitted. See [HI Theorem 3.2 p. 16] 
for instance and [U Theorem 8.3, p. 119] for a similar proof in the setting of upper 
triangular operators. We refer to [S3] for more on operator ranges and for their relations 
with de Branges Rovnyak spaces. 

7. Realizations of Schur multipliers 

The purpose of this section is to prove the following result, which is the counterpart 
in the quaternionic setting of the realization result for Schur functions. We first recall 
the following definition (see for instance [5], p. 14] in the complex case). 

Definition 7.1. Let M' he a right quaternionic Hilhert space. The pair {C,A) G 
Li{J^,IF) X L(^) is called closely outer connected if 

oo 

(7.1) y ran = J^. 

n=0 

The following is the counterpart of [3 Theorem 2.2.1, p. 49] in the setting of slice 
holomorphic functions. 

Theorem 7.2. Let s he a function from the open unit hall ofM into H. Then, s is a 
Schur multiplier if and only if there exist a right quaternionic Hilhert space Jif and a 
coisometric operator 

(7.2) (^^ : ^©H^^©e 
such that s can he written as a power series 

oo 
n=0 

where 

(7.3) = 
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Assume the pair [C, A) closely outer connected. Then, it is unique up to an isometry 
of right quaternionic Hilbert spaces. 

We note that the reahzation f l7.3p is called closely outer connected when the pair (C, A) 
is closely outer connected. 

Using formula fl2.12p we have for \\pA\\ < 1 

oo oo 

Y^P^CA^'^'B = = Ci.SR\p^\A)B, 

n=l n=l 

and so 

(7.4) s{p) = D + Ci.S^\p-\A)B = D +pC^{I -pA)-*B. 

The strategy of the proof is as follows. Let s be a Schur multiplier, and let J^{s) be 
the associated reproducing kernel quaternionic Hilbert space. As in the classical case, 
we want to show that Jif{s) is the state space, in an appropriate sense in the present 
setting, of a coisometric realization of s. We use the same method as in [5] (see in 
particular p. 50 there), suitably adapted to the non commutativity of H. One considers 
the reproducing kernel quaternionic right Hilbert space with reproducing kernel ks{p, q), 
and define the relation R (that is the right vector subspace of (^^(s) ©H) x (^^(s) ©H) 
defined as the right linear span of elements of the form 

f [ ks{p, q)qu\ f {ks{p,q)_^ h{p,0))u + h{p, 0)qv\ 1 
U /'V (s(g)-s(0))n + s(0)gt; J/' 

We claim that the relation R is densely defined and isometric. It will follow that R 
can be extended in a unique way to the graph of an isometric operator from J^{s) ©H 
into itself. This operator (or more precisely its adjoint) will give the realization. We 
first prove some preliminary lemmas and then give the proof of the theorem along the 
above lines. 

Lemma 7.3. The relation R is isometric and densely defined. 

Proof. We want to prove that 
(7.5) 

[ks{p, qi)qiui, ks{p, q2)q2U2].^(s) + [gi^i, q2V2]m = 
= [{ks{p, qi) - Kip, 0))mi + ksip, 0)qlvi, {ks{p, ^2) - ks{p, 0))u2 + ks{p, 0)q^V2].^(s)+ 
+ [(KgO - s(0))wi + 's(0)qlvi, - 7(0))U2 + '^q^V2]M 

for all choices of Ui,U2,Vi,V2 G and qi,q2 in the open unit ball of H. We rewrite 
this equality as 

u^MiUi + V2M2V1 + v^M^ui + u^MiVi = u^NiUi + V2N2V1 + v^N^^ui + u^N^vi, 

and we will show that Mj = Ni for i = 1,2,3,4. We write out in details the case i = 1, 
and only outline the other cases. 

Checking Mi = Nii Using the reproducing kernel property, we see that 

Ml = q2ks{q2,qi)qi- 
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Still with this same property, we have 

Ni = h{q2, qi) - h{q2, 0) - A;,(0, gi) + A;,(0, 0) + 



+ (.(g2)-3(0))(s(gi)-s(0)) 



= h{q2, qi) - (1 - s(g2)s(0)) - (1 - s(0)s(gi)) + 1 - \s{Q)Y 
+ sfe)^ - s{q2W) - s(0))^ + |s(0) p 
and thus the to prove the equality Mi = A'^i is equivalent to check that 



ks{,q2, qi) - q2ksiq2, qi)qi = l - s(g2)s(gi), 
but this is a direct consequence of the definition of the kernel kg. 

Checking M2 = N2: We now have M2 = q2qi and 

N2 = q2ks{0,0)q{ + q2\s{0)\% = q2(h = M2. 
We have M3 = M4 = and so we now need to check that N-^ = N4 = 0. 

Checking N3 = 0: This amounts to verify that 



q2{ks{0, qi) - ks{Q, 0)) + g2(s(0)(s(gi) - s(0)) = 0, 
but this is clear from the definition of the kernel since 

A:,(0,gi) = 1 - s(0)s(gi) and A:,(0, 0) = 1 - ' "^^^'^ 
Checking iV4 = 0: This amounts to verify that 



(A;,(0, 0) - A;,(0, ga))?! + (s(0) - s(g2))s(0)gr = 0, 
but this is also plain from the definition of the kernel kg- □ 
Lemma 7.4. R is the graph of a densely defined isometry. Let us denote by 

its extension to all of Jif{s) © H. Then 



(7.6) 



mip) = 


l/i, 


/(O)), 


p^O 
p = 0, 


{Bv){p) = 


\p-\s{p) - 
\si, 




p^O 
p = 0, 


Cf = 


/(O), 






Dv = 


s{0)v. 







Proof. We first check that R is the graph of a densely defined isometry. By definition, 
the domain of R is the set of F G ^(s) © EI such that there exists G G J^{s) © EI 
such that {F,G) G R. It is therefore dense by construction since the kernels ks{-,q)qu 
are dense in Jif{s). The isometry property implies that 

(0, G) G ^ G = 0. 
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We can thus introduce a densely defined operator T such that 

G = TF ^ (F, G) e R. 

T is a densely defined isometry since R is isometric. As in the case of complex Hilbert 
spaces, it extends to an everywhere defined isometry. We now compute the operator 
A. Let g G B and n G H. We have 

A*{k,{-, q)q)u = {ks{-, q) - ks{-, 0)) u. 
Hence, for / G J^{s) it holds that: 

m{Af){q) = [Af,ks{-,q)qu],^^s) 

= [f,ksi-,q) - ksi-,0)].jr-'{s) 
= n{f{q)-fm 

and hence 

q{Af){q) = f{q)-f{0). 

Similarly we have 

B*{h{;q)q)u = {J{q) -JiO)H 
so that we can write for v G H: 

uq{Bv){q) = [Bv, ks{-, q)u],%'(s) 

= [v, {s{q) - s(0))m]h 

= u{s{q) — s{0))v 

and hence the formula for B. To compute C we note that C*{qv) = ks{-, 0)qv for every 
g, f G H. So, for / G J^{s) we have: 

vqCf = [Cf,qv]m = [f, ksi-,0)qv],^(s) = vqf{0), 
and so Cf = /(O). Finally, it is clear that D = s(0). □ 
With these results at hand we turn to the proof of the realization theorem. 

Proof of Theorem 17. 2t We note that the pair (C, A) in (17. 6p is closely outer con- 
nected. Let / G J^{s), with power series 

oo 
n=0 

We have the formulas 

fn = CA-f, n = 0,l,2,... 

and 

f{p) = C^{I-pA)-y, fejr{s). 

Applying these formulas to the function 51 we obtain 

s{p) - s(0) =pCi.{I - pAy^Bl. 

We now turn to the converse and assume that s is given by a power series which 
converges in B and for which the coefficients are of the form (17. 3p . To prove that s is 
a Schur multiplier, we will check the formula 

(7.7) l-sip)J(q) = Uip)iUiq)r-pUip)iUiq)yq, p, g G B, 
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where 

oo 

(7.8) f/(p) = ^p"CA^ 

n=0 

We have 

oo oo 

1 - s{p)'^ =1-{D + J]p"CA"-i51)(D + J2 q'^CA'^-^Biy 

n=l m=l 

oo oo 

= 1 - DD* - DB*{A'^-^yC*f^ - 

?n=l n=l 

oo 
n,m=l 

oo oo 

= cc* + Y, c{A'yc*q"' + 

m=l n=l 

oo 

- J2 P'"CA''-\I ^ AA*){A"'-yc*q"' 

n,m=l 

= U{p){U{q)r-pU{p){U{q)rq, 

where U is as in fl7.8p and we have used the fact that the operator matrix f l7.2p is 
coisometric. 

It follows from (17. 7p that 

oo 

(7.9) 5^p"(l - s{p)-s{q))lt = U{p)U{q)\ 

n=0 

is positive definite in B, and so s is a Schur multiplier. 
Finally we turn to the uniqueness claim. Let 



and 

'A2 B2 
C2 D2 



be two closely outer-connected coisometric realizations of s, with state spaces right 
quaternionic Hilbert spaces Jifi and respectively. From (17. 9p we have 

Uiip)mq)y = u2ip)mp)r, 

where Ui and U2 are built as in (17.80 from the present realizations. It follows that 

CAUATYC; = C2Al{A^rc;, \/n,m e Nq. 
In view of the presumed outer-connectedness, the relation in Jifi x J^2 defined by 

{iAirC*u,iA;rc;u), uem, me No, 
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is a densely defined isometric relation with dense range. It is thus the graph of a 
unitary map U such that: 

U {{AirClu) = {A;rc;u, m e No, and n G H. 

Setting m = leads to t/CJ" = Q, that is 

(7.10) Ci = C2U. 
With this equality, writing 

(UADiiAirCi) = AliAirc; = A*UU*{A*rc; = {A*U){{AirCi), 

and taking into account that both pairs (Ci, Ai) and (C2, A2) are closely outer-connected, 
we obtain AiU* = U*A2, that is 

(7.11) UAi = A2U. 

Since clearly Di = D2 = s(0), it remains only to prove that UBi = B2. This follows 
from the equalities (where we use f l7.10p and (17.111) ) 

Sn = CiAI-^Bi = C2A^-^B2 = CiAl-^U*B2, n = 1, 2, . . . 

and from the fact that (Ci, Ai) is closely outer connected. □ 

We note that we have followed the arguments in [5] suitably adapted to the present 
case. In particular the proof of the uniqueness is adapted from that of [5l Theorem 
2.1.3, p. 46]. 

The preceding theorem gives a characterization of all Schur multipliers. A simple 
example is given by the choice 

(A B\ _ ( ^ vT^H^\ 

[C d)-[./T^ -a J' 

where a G B. The corresponding Schur multiplier Sa{p) is 

Sa{p) = -a + (1 - |ap)(l - pa)'* = {p - a) -k {1 - a)"*, 

and is the counterpart of an elementary Blaschke factor. The corresponding space 
J^{sa) is finite dimensional and spanned by the function (1 —pa)~*. We postpone to 
a future publication the study of finite dimensional (possibly indefinite) de Branges 
Rovnyak spaces in the present setting. 

We also note that (17.71) suggests an equivalent definition of Schur multiplier: The 
function s is a Schur multiplier if there is a function k{p, q) positive definite in B, and 
such that 

(7.12) 1 - s{p)s{q) = k{p,q) - pk{p,q)q, p, g G B, 

See also (19. ip . (19. 2p . These equations, as well as (I7.12p . can be rewritten as sums of 
positive definite functions, which induce an isometry. For instance, we can rewrite 

frm]) as 

1 + pk{p, q)q = k{p, q) + s{p)s{q), p,qeM, 
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We will not pursue this line of idea here, which is also called the lurking isometry 
method] see [22], and postpone it to a work where we consider functions of several 
quaternionic variables. 

Proposition 7.5. The reproducing kernel ks can be written as 

uh{p,q)v=[{C^{I-qAyyv, (C ^ (/ - 
(7.13) /~ \ 

\n=0 / 

Prooj. Indeed, both fc,(p,g) and U{p){U{q)y = (Er=oP"C'A"A*"C*g'^) satisfy the 
equation 



X — pXp = 1 — s{p)s{q). 

□ 



8. The Schur algorithm in the quaternionic setting 

Recall that we have denoted by S the set of functions analytic and contractive in the 
open unit disk D. We rewrite the recursion (11. 2p in a way which is used in signal 
processing, and is conducive to generalization to more general cases, and in particular 
to this case. We write 

(1 ?)(o ?) = (^(1 - ^(^)^) -(^W-^o)) 

^^•^^ = {1 - s{z)s^)-' {z -zs«(z)) 

= {i-s{z)-or'z{i -s«(^)). 

Let s{p) = J2'^=oP"'^ri be a Schur multiplier, and assume that |so| < 1. We now want 
the counterpart of (18. ip in the quaternionic setting. We set 

s{p) - s(0) =pa^^\p). 

Then, the counterpart of (18. ip in the quaternionic setting is: 

(/ -Ms)(^M_ ^")(^" °) = ((/ - M.M^)M, -(M.-M.J) 

= (/ - M.Mso) {Mp -{I - MsM^y^MpM^,) 
= {I- M,M^)Mp ^ (/ -(/ - M,iV%)-iM,J . 

The function s*^^^ defined by 

M,(i) = (/ - M,M^)-^M„, 

will be called the Schur transform of s. In view of Theorem 16.61 it is a Schur multiplier. 
When |s*-^''(0)| < 1, the preceding procedure can be iterated. This is the quaternionic 
counterpart of the Schur algorithm. 
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9. Concluding remarks 

Schur analysis has been extended to wide range of other settings. We mention in par- 
ticular the definitions of the Schur-Agler classes for functions analytic in the polydisk 
or the unit ball. See [H |2T1 [201 [221 [S2] • Other settings include the non- commutative 
case, [91123], Riemann compact surfaces [17] functions defined on an homogeneous tree, 
[18] , and others. See for instance [1] for the time- varying case, [10] for the stochastic 
case and [12] for the case of multi-dimensional systems. 

We mentioned in the paper a number of problems which will be considered in future 
publications. We here list some more: 

(1) We have presented a coisometric realization of a Schur multiplier. There are also iso- 
metric and unitary realizations, associated associated to the reproducing kernel Hilbert 
spaces ^^(5") and J^{Ds), where s" and Ds are defined in (15. 4p and (15. 5p respectively. 
Their counterpart in the setting of slice holomorphic functions will be presented else- 
where. 

(2) Let J e H"^" such that J = J* = J~\ A H"^"- valued function 6 defined in an 
open subset of B will be called J-contractive if the kernel 

00 

Keip, q) = J2p^{J - e(p) J(e(g))*)r = (1 - e(p) J(e(g))*) ^ (1 - pgr* 

n=0 

is positive definite. The study of these functions, and their relations to interpolation 
problems for Schur multipliers, are an important topic in Schur analysis. In the clas- 
sical case, these functions are called J-contractive, and are characteristic operator 
functions; see [551 EH ESI EDI EI]- Their structure has been described in [5U] . 

(3) We have mentioned Schur-Agler classes associated to the unit ball and the poly-disk. 
Their counterpart in the setting of Fueter series were presented in [161 113 • The study 
of Schur multipliers slice hyperholomorphic in several quaternionic variables should be 
a topic of interest. Let p = (pi, . . . , p^r) and q = (gi, . . . , qn) vary in B^. A Schur-Agler 
function is defined function s{p) such that the 

N 

(9.1) 1 - s{p)s{q) = ^{kj{p, q) - Pjkj{p, q)q]), (polydisk case), 
or 

N 

(9.2) 1 — s{p)s{q) = k{p, q) — '^^pjk{p, q)qj, (unit ball case), 
where the kj and k are positive definite kernels. 
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